Evolution and expression patterns of the neo-sex chromosomes of the crested ibis.

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Genomic evolution of island birds from the view of the Swinhoe's pheasant (Lophura swinhoii).

Island endemic birds account for the majority of extinct vertebrates in the past few centuries. To date, the evolutionary characteristics of island endemic bird's is poorly known. In this research, we de novo assembled a high-quality chromosome-level reference genome for the Swinhoe's pheasant, which is a typical endemic island bird. Results of collinearity tests suggest rapid ancient chromosome rearrangement that may have contributed to the initial species radiation within Phasianidae, and a role for the insertions of CR1 transposable elements in rearranging chromosomes in Phasianidae. During the evolution of the Swinhoe's pheasant, natural selection positively selected genes involved in fecundity and body size functions, at both the species and population levels, which reflect genetic variation associated with island adaptation. We further tested for variation in population genomic traits between the Swinhoe's pheasant and its phylogenetically closely related mainland relative the silver pheasant, and found higher levels of genetic drift and inbreeding in the Swinhoe's pheasant genome. Divergent demographic histories of insular and mainland bird species during the last glacial period may reflect the differing impact of insular and continental climates on the evolution of species. Our research interprets the natural history and population genetic characteristics of the insular endemic bird the Swinhoe's pheasant, at a genome-wide scale, provides a broader perspective on insular speciation, and adaptive evolution and contributes to the genetic conservation of island endemic birds.

A chromosome-level genome assembly of the Henosepilachna vigintioctomaculata provides insights into the evolution of ladybird beetles.

The ladybird beetle Henosepilachna vigintioctomaculata is an economically significant oligophagous pest that induces damage to many Solanaceae crops. An increasing number of studies have examined the population and phenotype diversity of ladybird beetles. However, few comparative genome analyses of ladybird beetle species have been conducted. Here, we obtained a high-quality chromosome-level genome assembly of H. vigintioctomaculata using various sequencing technologies, and the chromosome-level genome assembly was ~581.63 Mb, with 11 chromosomes successfully assembled. The phylogenetic analysis showed that H. vigintioctomaculata is a more ancient lineage than the other three sequenced ladybird beetles, Harmonia axyridis, Propylea japonica, and Coccinella septempunctata. We also compared positively selected genes (PSGs), transposable elements (TEs) ratios and insertion times, and key gene families associated with environmental adaptation among these ladybird beetles. The pattern of TEs evolution of H. vigintioctomaculata differs from the other three ladybird beetles. The PSGs were associated with ladybird beetles development. However, the key gene families associated with environmental adaptation in ladybird beetles varied. Overall, the high-quality draft genome sequence of H. vigintioctomaculata provides a useful resource for studies of beetle biology, especially for the invasive biology of ladybird beetles.

A chromosome-level genome of Semiothisa cinerearia provides insights into its genome evolution and control.

BACKGROUND: Semiothisa cinerearia belongs to Geometridae, which is one of the most species-rich families of lepidopteran insects. It is also one of the most economically significant pests of the Chinese scholar tree (Sophora japonica L.), which is an important urban greenbelt trees in China due to its high ornamental value. A genome assembly of S. cinerearia would facilitate study of the control and evolution of this species. RESULTS: We present a reference genome for S. cinerearia; the size of the genome was ~ 580.89 Mb, and it contained 31 chromosomes. Approximately 43.52% of the sequences in the genome were repeat sequences, and 21,377 protein-coding genes were predicted. Some important gene families involved in the detoxification of pesticides (P450) have expanded in S. cinerearia. Cytochrome P450 gene family members play key roles in mediating relationships between plants and insects, and they are important in plant secondary metabolite detoxification and host-plant selection. Using comparative analysis methods, we find positively selected gene, Sox15 and TipE, which may play important roles during the larval-pupal metamorphosis development of S. cinerearia. CONCLUSION: This assembly provides a new genomic resource that will aid future comparative genomic studies of Geometridae species and facilitate future evolutionary studies on the S. cinerearia.

Genome-wide survey of open chromatin regions in two swallowtail butterflies Papilio machaon and P. bianor.

Papilio machaon was assigned as the type species for all butterflies by Linnaeus and P. bianor is a congener but exhibits a great difference in morphology (especially larva and adult color pattern) and larval host plants from P. machaon. Thus, they are the ideal models to investigate genetic mechanisms underlying morphology and plasticity between congeners. The reference genomes of both species were dissected in our previous studies, but little is known about their regulatory genome and the epigenetic regulation of gene expression throughout developmental stages. Here, we profiled the chromatin accessibility and gene expression of three developmental stages (the 4th instar larva [L4], the 5th instar larva [L5], and pupa [P]) using transposase accessible chromatin sequencing (ATAC-seq) and RNA-seq. Results showed that many accessible chromatin peaks were identified at three developmental stages (peak number, P. machaon: 44,977 [L4], 36,919 [L5], 47,147 [P]; P. bianor: 20,341 [L4], 44,668 [L5], 62,249 [P]). Moreover, the number of differentially accessible peaks and differentially expressed genes between larval stages of each butterfly species are significantly fewer than that between larval and pupal stages, suggesting a higher similarity within larvae and a significant difference between larvae and pupae. This study added the annotated information of chromatin accessibility genome-wide of the two papilionid species and will promote the investigation of gene regulation in butterfly evolution.

A chromosome-level genome of Brachymystax tsinlingensis provides resources and insights into salmonids evolution.

Brachymystax tsinlingensis Li, 1966 is an endangered freshwater fish with economic, ecological, and scientific values. Study of the genome of B. tsinlingensis might be particularly insightful given that this is the only Brachymystax species with genome. We present a high-quality chromosome-level genome assembly and protein-coding gene annotation for B. tsinlingensis with Illumina short reads, Nanopore long reads, Hi-C sequencing reads, and RNA-seq reads from 5 tissues/organs. The final chromosome-level genome size is 2,031,709,341 bp with 40 chromosomes. We found that the salmonids have a unique GC content and codon usage, have a slower evolutionary rate, and possess specific positively selected genes. We also confirmed the salmonids have undergone a whole-genome duplication event and a burst of transposon-mediated repeat expansion, and lost HoxAbß Hox cluster, highly expressed genes in muscle may partially explain the migratory habits of B. tsinlingensis. The high-quality B. tsinlingensis assembled genome could provide a valuable reference for the study of other salmonids as well as aid the conservation of this endangered species.

Chromosome-level genome assembly of the black widow spider Latrodectus elegans illuminates composition and evolution of venom and silk proteins.

BACKGROUND: The black widow spider has both extraordinarily neurotoxic venom and three-dimensional cobwebs composed of diverse types of silk. However, a high-quality reference genome for the black widow spider was still unavailable, which hindered deep understanding and application of the valuable biomass. FINDINGS: We assembled the Latrodectus elegans genome, including a genome size of 1.57 Gb with contig N50 of 4.34 Mb and scaffold N50 of 114.31 Mb. Hi-C scaffolding assigned 98.08% of the genome to 14 pseudo-chromosomes, and with BUSCO, completeness analysis revealed that 98.4% of the core eukaryotic genes were completely present in this genome. Annotation of this genome identified that repetitive sequences account for 506.09 Mb (32.30%) and 20,167 protein-coding genes, and specifically, we identified 55 toxin genes and 26 spidroins and provide preliminary analysis of their composition and evolution. CONCLUSIONS: We present the first chromosome-level genome assembly of a black widow spider and provide substantial toxin and spidroin gene resources. These high-qualified genomic data add valuable resources from a representative spider group and contribute to deep exploration of spider genome evolution, especially in terms of the important issues on the diversification of venom and web-weaving pattern. The sequence data are also firsthand templates for further application of the spider biomass.

Genomic deciphering of sex determination and unique immune system of a potential model species rare minnow (Gobiocypris rarus).

Gobiocypris rarus is sensitive to environmental pollution, especially to heavy metal and grass carp reovirus (GCRV). Hence, it has potential utility as a biological monitor. Genetic deciphering of its unique immune system will advance our understanding of its unique adaptive strategies, which provide cues for its better application. A de novo genome of rare minnow was obtained, and its sex determination mechanism is ZZ/ZW. We identified several specific mutation genes and specific lost genes of rare minnow, and these might be related to the sensitivity of rare minnow to environmental stimuli. We also analyzed the gene expression level of different organs/tissues and found that several IFIT genes may play key roles in GCRV resistance. In addition, knockout of the gene PCDH10L indicates that PCDH10L affects Pb2+-induced mortality in rare minnow. Rare minnow is ready for genetic manipulation and shows potential as an emerging experimental model.

Genomic insights into the evolution of the critically endangered soft-shelled turtle Rafetus swinhoei.

The 2019 death of the female Swinhoes soft-shelled turtle (Rafetus swinhoei) individual in China put this species in a dire conservation state. To avoid the loss of genetic information from the female individual and facilitate the protection and preservation of R. swinhoei, we sequenced and assembled its draft chromosome-level genome by combining BGISEQ short reads, Nanopore long reads, and Hi-C sequencing reads. We characterized its potential ZZ/ZW sex determination system. Phylogeny and divergence time analyses indicated that R. swinhoei has a close relationship with Pelodiscus sinensis, and they shared a common ancestor ~54.4 Ma. Population history analyses indicated that R. swinhoei has the lowest heterozygosity among all turtles examined, and its population size has declined dramatically over the past 10 million years. Many core genes involved in autophagy and DNA damage response, such as Rab5a and Parp1, were under positive selection and/or evolving rapidly, and these genes may make important contributions to the long lifespan and evolutionary adaptation of R. swinhoei. The loss of several tooth-related genes in turtle genomes explains the genetic basis of their toothless phenotype. This study not only provides a genomic resource for scientists to study the species using the strategy of comparative genomics, but also provides a data basis for us to find more potential R. swinhoei individuals in the wild.

Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms.

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

Chromatin accessibility profiling provides insights into larval cuticle color and adult longevity in butterflies.

Butterflies are diverse in virtually all aspects of their ontogeny, including morphology, life history, and behavior. However, the developmental regulatory mechanisms underlying the important phenotypic traits of butterflies at different developmental stages remain unknown. Here, we investigated the developmental regulatory profiles of butterflies based on transposase accessible chromatin sequencing (ATAC-seq) at three developmental stages in two representative species ( Papilio xuthus and Kallima inachus). Results indicated that 15%-47% of open chromatin peaks appeared in associated genes located 3 kb upstream (i.e., promoter region) of their transcription start site (TSS). Comparative analysis of the different developmental stages indicated that chromatin accessibility is a dynamic process and associated genes with differentially accessible (DA) peaks show functions corresponding to their phenotypic traits. Interestingly, the black color pattern in P. xuthus 4th instar larvae may be attributed to promoter peak-related genes involved in the melanogenesis pathway. Furthermore, many longevity genes in 5th instar larvae and pupae showed open peaks 3 kb upstream of their TSS, which may contribute to the overwintering diapause observed in K. inachus adults. Combined with RNA-seq analysis, our data demonstrated that several genes enriched in the melanogenesis and longevity pathways also exhibit higher expression, confirming that the expression of genes may be closely related to their phenotypic traits. This study offers new insights into larval cuticle color and adult longevity in butterflies and provides a resource for investigating the developmental regulatory mechanisms underlying butterfly ontogeny.

Large-scale sequencing of flatfish genomes provides insights into the polyphyletic origin of their specialized body plan.

The evolutionary and genetic origins of the specialized body plan of flatfish are largely unclear. We analyzed the genomes of 11 flatfish species representing 9 of the 14 Pleuronectiforme families and conclude that Pleuronectoidei and Psettodoidei do not form a monophyletic group, suggesting independent origins from different percoid ancestors. Genomic and transcriptomic data indicate that genes related to WNT and retinoic acid pathways, hampered musculature and reduced lipids might have functioned in the evolution of the specialized body plan of Pleuronectoidei. Evolution of Psettodoidei involved similar but not identical genes. Our work provides valuable resources and insights for understanding the genetic origins of the unusual body plan of flatfishes.

The Genomes of Two Billfishes Provide Insights into the Evolution of Endothermy in Teleosts.

Endothermy is a typical convergent phenomenon which has evolved independently at least eight times in vertebrates, and is of significant advantage to organisms in extending their niches. However, how vertebrates other than mammals or birds, especially teleosts, achieve endothermy has not previously been fully understood. In this study, we sequenced the genomes of two billfishes (swordfish and sailfish), members of a representative lineage of endothermic teleosts. Convergent amino acid replacements were observed in proteins related to heat production and the visual system in two endothermic teleost lineages, billfishes and tunas. The billfish-specific genetic innovations were found to be associated with heat exchange, thermoregulation, and the specialized morphology, including elongated bill, enlarged dorsal fin in sailfish and loss of the pelvic fin in swordfish.

African lungfish genome sheds light on the vertebrate water-to-land transition.

Lungfishes are the closest extant relatives of tetrapods and preserve ancestral traits linked with the water-to-land transition. However, their huge genome sizes have hindered understanding of this key transition in evolution. Here, we report a 40-Gb chromosome-level assembly of the African lungfish (Protopterus annectens) genome, which is the largest genome assembly ever reported and has a contig and chromosome N50 of 1.60 Mb and 2.81 Gb, respectively. The large size of the lungfish genome is due mainly to retrotransposons. Genes with ultra-long length show similar expression levels to other genes, indicating that lungfishes have evolved high transcription efficacy to keep gene expression balanced. Together with transcriptome and experimental data, we identified potential genes and regulatory elements related to such terrestrial adaptation traits as pulmonary surfactant, anxiolytic ability, pentadactyl limbs, and pharyngeal remodeling. Our results provide insights and key resources for understanding the evolutionary pathway leading from fishes to humans.

Genome-wide analysis of DNA G-quadruplex motifs across 37 species provides insights into G4 evolution.

G-quadruplex (G4) structures have been predicted in the genomes of many organisms and proven to play regulatory roles in diverse cellular activities. However, there is little information on the evolutionary history and distribution characteristics of G4s. Here, whole-genome characteristics of potential G4s were studied in 37 evolutionarily representative species. During evolution, the number, length, and density of G4s generally increased. Immunofluorescence in seven species confirmed G4s' presence and evolutionary pattern. G4s tended to cluster in chromosomes and were enriched in genetic regions. Short-loop G4s were conserved in most species, while loop-length diversity also existed, especially in mammals. The proportion of G4-bearing genes and orthologue genes, which appeared to be increasingly enriched in transcription factors, gradually increased. The antagonistic relationship between G4s and DNA methylation sites was detected. These findings imply that organisms may have evolutionarily developed G4 into a novel reversible and elaborate transcriptional regulatory mechanism benefiting multiple physiological activities of higher organisms.

Chromosome-level genome assembly of Lethenteron reissneri provides insights into lamprey evolution.

The reissner lamprey Lethenteron reissneri, belonging to the class Cyclostomata, serves as a bridge between invertebrates and jawed vertebrates, and is considered the sister group of jawed vertebrates. However, despite this evolutionary significance, the genetic mechanisms underlying the adaptive evolution of the lamprey lineage remain unclear. Here, we assembled a 1.06 Gb chromosome-level draft genome of L. reissneri, with 72 chromosomes (ranging in length from 4.5 Mb to 25.9 Mb) and a scaffold N50 length of 13.23 Mb. Genome quality comparisons revealed that the reissner lamprey genome has higher completeness and contiguity than the previously published sea lamprey and Japanese lamprey genomes. Moreover, reissner lamprey, sea lamprey, and Japanese lamprey species share similar transposable element profiles and Hox gene cluster compositions, suggesting that a burst of transposable element activity and whole genome duplication occurred before their divergence. Additionally, the Lip gene copy numbers, which have been studied for their functions in the host defence system, were found to be expanded uniquely in lamprey lineages, suggesting key roles for these genes in lamprey evolution and adaptation. We also identified two neural-related genes, Nrn1 and Unc13a, with copy number expansions in jawed vertebrates, which may be functionally relevant to the origin of lamprey brains. Hence, this study not only provides the first chromosome-level reference genome for Cyclostomata, but also highlights features of the unique biology and adaptive evolution of the lamprey lineage.

Chromosome-level genome assembly of Paralithodes platypus provides insights into evolution and adaptation of king crabs.

The blue king crab, Paralithodes platypus, which belongs to the family Lithodidae, is a commercially and ecologically important species. However, a high-quality reference genome for the king crab has not yet been reported. Here, we assembled the first chromosome-level blue king crab genome, which contains 104 chromosomes and an N50 length of 51.15 Mb. Furthermore, we determined that the large genome size can be attributed to the insertion of long interspersed nuclear elements and long tandem repeats. Genome assembly assessment showed that 96.54% of the assembled transcripts could be aligned to the assembled genome. Phylogenetic analysis showed the blue king crab to have a close relationship with the Eubrachyura crabs, from which it diverged 272.5 million years ago. Population history analyses indicated that the effective population of the blue king crab declined sharply and then gradually increased from the Cretaceous and Neogene periods, respectively. Furthermore, gene families related to developmental pathways, steroid and thyroid hormone synthesis, and inflammatory regulation were expanded in the genome, suggesting that these genes contributed substantially to the environmental adaptation and unique body plan evolution of the blue king crab. The high-quality reference genome reported here provides a solid molecular basis for further study of the blue king crab's development and environmental adaptation.

Maintenance treatment of transcranial magnetic stimulation (TMS) for treatment-resistant depression patients responding to acute TMS treatment.

A growing body of studies has demonstrated that acute transcranial magnetic stimulation (TMS) therapy for treatment-resistant major depressive disorder (MDD) has achieved significant antidepressant effects and can alleviate other related symptoms. However, MDD has a high relapse rate, and patients with depressive symptoms can relapse weeks or months after acute TMS treatment. The lack of necessary TMS maintenance protocols after completing acute TMS treatment with full remission might be one of the reasons for the high relapse rates in MDD patients. Thus, investigating post-TMS treatment maintenance guidelines is important for decreasing relapse in treatment-resistant depression patients who had initially responded to acute TMS therapy. Therefore, we recommend a scientific approach to decrease relapse in treatment-resistant depression patients who had initially responded to acute TMS treatment.

Genomic and experimental data provide new insights into luciferin biosynthesis and bioluminescence evolution in fireflies.

Fireflies are among the most charismatic insects for their spectacular bioluminescence, but the origin and evolution of bioluminescence remain elusive. Especially, the genic basis of luciferin (D-luciferin) biosynthesis and light patterns is largely unknown. Here, we present the high-quality reference genomes of two fireflies Lamprigera yunnana (1053 Mb) and Abscondita terminalis (501 Mb) with great differences in both morphology and luminous behavior. We sequenced the transcriptomes and proteomes of luminous organs of two species. We created the CRISPR/Cas9-induced mutants of Abdominal B gene without luminous organs in the larvae of A. terminalis and sequenced the transcriptomes of mutants and wild-types. Combining gene expression analyses with comparative genomics, we propose a more complete luciferin synthesis pathway, and confirm the convergent evolution of bioluminescence in insects. Using experiments, the function of the firefly acyl-CoA thioesterase (ACOT1) to convert L-luciferin to D-luciferin was validated for the first time. Comparisons of three-dimension reconstruction of luminous organs and their differentially expressed genes among two species suggest that two positive genes in the calcium signaling pathway and structural difference of luminous organs may play an important role in the evolution of flash pattern. Altogether, our results provide important resources for further exploring bioluminescence in insects.